Casimir effect due to a single boundary as a manifestation of the Weyl problem 
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The Casimir self-energy of a boundary is ultraviolet-divergent. In many cases the divergences can 
be eliminated by methods such as zeta-function regularization or through physical arguments (ultra- 
violet transparency of the boundary would provide a cutoff) . Using the example of a massless scalar 
field theory with a single Dirichlet boundary we explore the relationship between such approaches, 
with the goal of better understanding the origin of the divergences. We are guided by the insight due 
to Dowker and Kennedy (1978) and Deutsch and Candelas (1979), that the divergences represent 
measurable effects that can be interpreted with the aid of the theory of the asymptotic distribution 
of eigenvalues of the Laplacian discussed by Weyl. In many cases the Casimir self-energy is the sum 
of cutoff-dependent (Weyl) terms having geometrical origin, and an "intrinsic" term that is inde- 
pendent of the cutoff. The Weyl terms make a measurable contribution to the physical situation 
even when regularization methods succeed in isolating the intrinsic part. Regularization methods 
fail when the Weyl terms and intrinsic parts of the Casimir effect cannot be clearly separated. 
Specifically, we demonstrate that the Casimir self-energy of a smooth boundary in two dimensions 
is a sum of two Weyl terms (exhibiting quadratic and logarithmic cutoff dependence), a geometrical 
term that is independent of cutoff, and a non-geometrical intrinsic term. As by-products we resolve 
the puzzle of the divergent Casimir force on a ring and correct the sign of the coefficient of linear 
tension of the Dirichlet line predicted in earlier treatments. 

PACS numbers: 03.70.+k, ll.10.-z, ll.10.Gh, 42.50.Pq 
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I. INTRODUCTION 

Casimir interactions are the macroscopic response of 
the physical vacuum to the introduction of boundaries. 
They were first derived as an attractive force between 
perfectly conductive parallel plates induced by the zero- 
point motion of the electromagnetic field There is 
convincing experimental evidence for the reality of these 
forces [2| and a vast body of literature dedicated to var- 
ious aspects of the phenomenon Q. 

The Casimir interaction £ is the difference between the 
vacuum energy of the system constrained by the bound- 
aries and that of free space. Since boundaries made of 
real materials are transparent to sufficiently high-energy 
modes, the high-energy spectrum is unaffected by the 
geometry of the system, and only a finite range of the 
spectrum needs be considered [H|j|- However, in the the- 
oretical treatments of this effect the vacuum energies are 
usually calculated from an effective low-energy harmonic 
field theory (such as quantum electrodynamics in the case 
of the electromagnetic Casimir effect), so that they are 
approximated by the sum of zero-point energies of a col- 
lection of simple harmonic oscillators with a spectrum 
lj = c|k| (where c is the speed of light). In this model, 
the dispersion relation holds for arbitrarily large wave 
vectors k; both the "constrained" and "free" vacuum en- 
ergy densities are ultraviolet divergent, and the Casimir 
interaction is the difference between two infinite quanti- 
ties. This problem can be resolved by a soft-cutoff mod- 
ification of the large-k part of the spectrum that leads 
to a finite vacuum energy. When the result is not sensi- 
tive to the form of the cutoff, the Casimir interaction can 
be extracted by taking the cutoff to infinity at the end 



of calculation. However, it is more common that there 
are contributions that diverge in the absence of a cut- 
off; we will refer to these as "formally divergent terms," 
since on physical grounds the Casimir energy must be 
finite. Some approaches to the calculation make use of 
analytic continuation Q or dimensional regularization Q 
to eliminate the formally divergent terms and extract the 
"intrinsic" part that is independent of the cutoff. 

The important virtue of the model is that for many 
geometries the various approaches give the same result 
for the intrinsic part of the Casimir energy, that is in- 
dependent of the cutoff and depends only on Planck's 
constant H, the speed of light c, and macroscopic length 
scales. However, this is not always the case: in particular 
the divergences occurring for spherical geometry in even 
space dimensions are not removed by any regularization 
scheme @, H[ ■ What this means physically represents an 
open problem; it seems to imply that in two dimensions 
a conducting ring placed in vacuum is subject to an infi- 
nite stress. The two-dimensional version of this problem 
(a circle) was first studied by Sen Q , who concluded that 
in addition to the expected cutoff-independent terms, the 
Casimir energy contains contributions that depend on the 
cutoff. He further found that the formally divergent con- 
tributions can be interpreted geometrically. 

The main goal of this paper is to analyze the formal 
divergences encountered in the Casimir self-energy of an 
object with curved boundary. This is done by looking at 
the case of a massless scalar field theory with a Dirichlet 
boundary and explicitly taking into consideration the in- 
visibility of the boundary to sufficiently high-energy field 
modes. Specifically we a re g uided by observation due 
to Dowker and Kennedy [lj| and Deutsch and Cande- 
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las [IT| that the formally divergent terms correspond to 
measurable effects, and that they can be understood and 
interpreted geometrically through the connection of the 
Casimir problem to the Weyl problem of the asymptotic 
distribution of the eigenvalues of the Laplacian [12]. This 
connection will be illustrated for several geometries by si- 
multaneous solution of the Casimir and Weyl problems 
and showing that the Weyl piece of the Casimir effect 
corresponds to the formally divergent contributions. The 
cases for which it has not been possible to regularize the 
divergence are shown to belong to a separate class where 
the Weyl piece of the Casimir energy contains a contribu- 
tion that is logarithmically dependent on the cutoff. We 
argue that this is what happens in even space dimensions 
and demonstrate our contention by detailed analysis of 
the Casimir effect due to a smooth Dirichlet boundary in 
two spatial dimensions. As by-products we resolve the 
problem of the divergent Casimir force exerted on the 
ring in two dimensions 0, H[ and correct the sign of the 
coefficient of linear tension given previously [9] . 

The organization of this paper is as follows. In Section 
II we introduce a new computational method in which the 
Casimir effect is viewed as due to the vacuum fluctuations 
eliminated by introduction of the Dirichlet boundaries. 

In Section III we review the relationship between the 
Weyl problem of asymptotic distribution of the eigen- 
values of the Laplacian and the Casimir problem. As 
a motivation we calculate the formally divergent Casimir 
energy of a single Dirichlet plane (Section IIIA) and show 
how the latter can be understood geometrically in terms 
of the areal density of states (DOS) of the Weyl prob- 
lem (Section IIIB). For a field confined to a region the 
Weyl DOS is defined as a smooth part of the exact DOS. 
Since the Weyl DOS contains ultraviolet spectral infor- 
mation, it is natural to conjecture that the Weyl DOS 
is entirely responsible for the formally divergent parts of 
the Casimir effect. This leads to a large wave vector ex- 
pansion of the Weyl DOS. The Weyl-Casimir correspon- 
dence allows us to deduce the two leading terms of this 
expansion known as Weyl's conjecture. These ideas are 
illustrated in Section IIIC where we compute the Casimir 
energies of a field confined to a one-dimensional interval 
with Dirichlet and periodic boundary conditions. Ad- 
ditionally we demonstrate that the Weyl piece of the 
Casimir effect is insensitive to the topology of the prob- 
lem (whereas the intrinsic part is sensitive to this). The 
analysis of the same problem via the zeta-function reg- 
ularization technique carried out in Section HID reveals 
that the Weyl piece of the effect is usually discarded (but 
possibly can be retrieved). One of the consequences of 
the separability of the Weyl and intrinsic parts of the 
Casimir effect discussed in Section HIE is an expression 
for the Weyl energy of the boundary as a surface integral 
of an even function of curvature invariants. For the spe- 
cial case of a spherical shell in general space dimension, 
this expression becomes a finite series in powers of the 
sphere radius. Dimensional analysis of this series reveals 
the special role of even space dimensions as being cases 



where there is not a clear separation between the Weyl 
and cutoff-independent pieces of the Casimir energy. 

In Section IV we compute the Casimir energy in two 
spatial dimensions due to an arbitrary smooth Dirichlet 
boundary; this contains the main results of our work. 
This is the case where regularization techniques fail to 
give a finite value for the Casimir energy. As an illus- 
tration of this statement we show why the zeta-function 
regularization approach cannot be successfully employed 
here. Consistent with the observations made in Section 
HIE, we demonstrate that even though the Weyl-Casimir 
correspondence continues to play a prominent role, we no 
longer have a clear cut separation of the cutoff-dependent 
and intrinsic parts of the Casimir effect. Our result is an 
expression for the Casimir energy due to a smooth bound- 
ary, in the form of a sum of formally divergent terms com- 
ing from the Weyl expansion, a geometrical term that is 
cutoff independent, and a non-geometrical intrinsic con- 
tribution. Apart from the last term that cannot be given 
in general form, the bulk of the effect has the form of 
a purely geometrical expression. The latter is deduced 
from analysis of a particular case (a wedge bounded by 
an arc) for the Dirichlet and periodic boundary condi- 
tions (Sections IVA and IVB, respectively) that can be 
generalized to the case of an arbitrary curve. In Sec- 
tion IVC these results are applied to compute the cutoff- 
independent difference in the Casimir energies between 
the Dirichlet and periodic geometries. 

Our conclusions are summarized in Section V. For all 
the cases considered in this work we correctly predict the 
sign of the effect, which is possible because we calculate 
the Casimir energy in terms of the field modes elimi- 
nated by the Dirichlet boundary. Prediction of the sign 
of the Casimir effect is an important fundamental prob- 
lem whose general solution is yet to be found. The only 
rigorous result in this area is a proof [l3[ that the Casimir 
force between two arbitrarily shaped bodies related by re- 
flection is always attractive. Additionally, in the scalar 
case the sign of the effect can be estimated from the sign 
of the contribution due to shortest periodic rays 



II. CALCULATING THE CASIMIR ENERGY 

Here we will explain how to calculate the change in en- 
ergy resulting from the introduction of a Dirichlet bound- 
ary into a quantum system at zero temperature. For 
technical reasons, it is useful to first formulate the prob- 
lem at finite temperature T, and subsequently take the 
limit T — > 0. The starting place is a Gaussian field theory 
having the Euclidean action 

Se[w] = \J^ drd d x (V 2 (§^) 2 + (V W ) 2 ) . (1) 

The real scalar w is a function of the <i-dimcnsional posi- 
tion vector r and is periodic in imaginary time r, so that 
w(r,0) = w(r,h/T) [3. 
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The scalar field theory (p} can be viewed as a toy ver- 
sion of electricity and magnetism. The divergences en- 
countered for curved geometries exist in both theories, so 
they are not due to specifics of electricity and magnetism. 
Therefore we study the problem in the simpler setting of 
the scalar theory. 

The zero-point energy can be calculated by means of 
a functional integral that makes use of the correspon- 
dence between the Feynman path integral for the d- 
dimensional field theory and the partition function for 
a d + 1-dimensional classical statistical mechanics prob- 
lem. The path integral is 

Z w = J Dw(T,T)exp{-S E [w]/h). (2) 

This resembles the partition function for a classical sta- 
tistical mechanics problem, with the Hamiltonian re- 
placed by the action Se and the temperature replaced 
by Planck's constant h (please note that this is distinct 
from the real temperature T) [l6j]. In this analogy the 
energy for the quantum system corresponds to the "free 
energy" per unit "length" in the imaginary time direc- 
tion, so that £ = h(\nZ w )/(h/T) = -TlnZ w , and the 
zero-point energy is the T — > limit of this. 

Now assume that the vacuum is disturbed by the ad- 
dition of sharp boundaries which constrain the field in 
some way and thus eliminate some degrees of freedom of 
the vacuum fluctuations. The constrained field (which 
we will refer to as v) inherits any boundary conditions 
imposed on w as well as new conditions on surfaces Di 
of Dirichlet (v\u i = 0) type where the subscript i la- 
bels the boundaries. We will write the difference be- 
tween the original and constrained fields in the form 
w(r, t) = v(r, t) + u(r, r), where u(r, r) can be chosen to 
satisfy the d + 1-dimensional equation 

(-A_+A)t t = 0, u\ Dt =fi(T,T) (3) 

where fi are functions defined on the boundaries and 
determined by the boundary values of w; they play the 
role of dynamical variables of our approach. The reason 
for defining u this way is that it eliminates the cross term 
in the action, so that S^w] = S E [v] + Se[u}. Then the 
functional integral factors into integrations over v and u, 
so that Z w — Z V Z U . Since the zero-point energy for the 
unconstrained system is determined by Z w and the zero- 
point energy for the constrained system is determined by 
Z v , the Casimir energy is given by T — > limit of £ = 
const + Tin Z u . This result allows us to calculate directly 
the change in energy due to the fluctuation modes that 
have been eliminated. 

A simplification is achieved by expanding all the dy- 
namical variables of the problem into a Fourier series 
in the imaginary time domain; for example u(r, r) = 
tto)( r ) exp iujt where the Fourier coefficients u u (r) 
are solutions to the boundary-value problem for the 



Helmholtz equation 

(A-— K = 0, u u \ Di =/ w ,i(r) (4) 
c l 

The calculation of the action Se[u] is further simpli- 
fied when the identity (Vu) 2 = div(uS/u) — uAu is sub- 
stituted into Eq.([T]). Then the integral of div(uVu) over 
d d x transforms into a sum of surface integrals. The re- 
maining integral over dr vanishes due to the relation 
Au = —d 2 u/c 2 dr 2 and the condition of periodicity, 
u(r, 0) = u(r, h/T). As a result we find 

S E [u] = - J dT z2j [wVu] l ds l 

= 2TAZ fuA^U-ukdSi (5) 

ui,i 

Here [ip]i stands for the discontinuity of ip across the 
i-th boundary, the summation is performed over all the 
boundaries and we employed continuity of the field u = f 
at the Dirichlet boundary. 

Although our approach is applicable to an arbitrary 
number of the Dirichlet boundaries, in all the cases con- 
sidered in this paper only one boundary contributes into 
the integral ([S]). Since we are dealing with a harmonic 
field theory, the solution to the boundary-value problem 
(0| is linear in the surface field f u and thus the action 
Se is a quadratic diagonal form of f u : 

S e = 7^ ^ x 7TT7 ^ ( 6 ) 
2T ' A„( \uj \ c) 

where the subscript v (representing one or more indices) 
labels the normal modes of the field u that have been 
eliminated by the boundary in question, and \ v (\ui\/c) > 
is a set of geometry-dependent coefficients. Then ac- 
cording to the prescription outlined above, the Casimir 
energy is given by 

£ = const + -^ln(27rTA !/ (|w|/c)) 

-> const H V ' / du\n{2nT\ v {io/c)) (7) 

where in the second step we took the zero-temperature 
limit according to the rule — > (h/T) J du/2-K. As 
was already mentioned, any physical boundary is invis- 
ible to high-energy modes and from now on this fact is 
made implicit by adding a prime to all pertinent sums 
and integrals over normal modes. This means that large 
energy contributions into the sum in ([7]) are suppressed 
by a smooth cutoff function that represents transmission 
properties of the boundary. We will not need an explicit 
form of the cutoff function. 

The constant in (J7J is such as it cancels the temper- 
ature dependence and makes the argument of the log- 
arithm dimensionless. Rigorous determination of this 
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constant requires an analysis of the functional integra- 
tion measure in © that led to ([7} • Instead we provide a 
heuristic argument that gives the same answer. Since in 
our approach the Casimir effect arises due to harmonic 
field oscillators eliminated by the Dirichlet boundary it 
would suffice to understand the case of one degree of free- 
dom (the harmonic oscillator) and then the result can be 
generalized to that of the harmonic field theory (JTJ) . 

Let us consider a harmonic oscillator of mass m and 
natural frequency fi whose position is given by the dis- 
placement function x(r). The Euclidian action 



S B [x] = 'dr((g) 2 + nV 



mft 
2T 



J2(^ 2 + n 2 )\ x {uj)\ 



(8) 



is Gaussian and the zero-point energy is given by 

£ = const+— dcolnl 9 
2vr J Q \2TiT/mJ 

We need to determine the additive constant, and also deal 
with the divergence of the integral at the upper limit. 
Since the physical source of the zero-point energy is the 
confinement of the field, Eq. © should vanish in the free 
particle limit w/f2 — > oo. Then subtracting from Eq.© 
its high-frequency limit we obtain 
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2^ 



duj In 



n 2 



Ml 
~2~ 



(10) 



as it should be. A rigorous argument based on care- 
ful analysis of the functional integration measure leading 
from © to |10| has been given by Kleinert [17| . 

In the Casimir problem the confined field is represented 
by a collection of harmonic oscillators, and the same rule 
applied to Eq.© gives our final general expression for 
the Casimir energy: 



diu In 



duj hcu d ^ A„(cj/c) 
7r 2 duj A„(oo) 



(11) 



The infinity in the argument of \ v is not to be taken 
literally: in any problem the variable u would enter via 
the dimensionless combination uj/c times a macroscopic 
length scale. The infinity means that the limit of this 
length scale going to infinity is to be taken so that the 
zero-point energy is the confinement energy relative to 
a reference free field geometry. Our prescription (| 1 1|) 
agrees with that of Brevik and Elizalde [H] employed 
in their calculation of the Casimir energy of a piecewise 
uniform string. 

The rule (fTTj) guarantees that the integral over u> is 
ultraviolet convergent and that every field oscillator of 
frequency f2„ eliminated by the boundary contributes 



— ftf2„/2 into the Casimir energy. The second represen- 
tation of (TTTj) (where we integrated by parts) hints at a 
relationship between our approach and the calculation of 
the Casimir energy by means of a contour integral . 

We note that the integration variable uj in Eq. (JTTJ) has 
its origin in the Fourier transform in the imaginary time 
domain that led from Eq.© to Eq.©. This auxiliary 
variable is not related to the physical frequency of a field 
oscillator and therefore it should not be subject to any 
large-w cutoff. On the other hand, the mode index v 
is a physical variable akin to a wave vector and thus it 
should be subject to a cutoff as a consequence of the 
penetrability of the boundary to high-energy modes. 

Our cutoff rule is different from that employed by Sen 
in his analysis of the circle geometry. Expressed in 
our language, he imposes an exponential cutoff on the 
variable oj but leaves the mode variable v unaffected. Al- 
though his procedure gives a finite answer for the circular 
geometry, this would not be the case for the cylinder and 
many other geometries. On the other hand, our rule (|11[) 
is universally applicable. 



III. GEOMETRICAL INTERPRETATION OF 
ULTRAVIOLET DIVERGENCES 

Let us assume that the physical boundary is charac- 
terized by a frequency cutoff at scale luq. The boundary 
is impenetrable to low-energy field modes but invisible 
to the modes of energy significantly greater than ftwo- 
Such a boundary can be modeled by a Dirichlet surface. 
We will determine the coefficient 70 of the fluctuation- 
induced surface tension of a single Dirichlet plane im- 
mersed in a d-dimensional vacuum. Since this problem 
is only characterized by microscopic energy and length 
scales, ftwo and c/loq, respectively, the outcome can be es- 
timated via dimensional analysis. Indeed the only quan- 
tity having dimension of energy divided by length to the 
power d — 1 is 70 ~ hu a / {c/loq)^ 1 = hc(uj Q /c) d . We 
observe that the surface tension diverges as the cutoff 
frequency coq approaches infinity. This divergence has no 
effect on the Casimir stress on parallel planes because 
the overall area of the boundaries remains fixed as the 
distance between the planes is varied. 

However this is not the case for curved boundaries. 
For example, a change of the radius of a circle implies a 
change of the perimeter and as a result the formally diver- 
gent self-energy will contribute to the Casimir force per 
length. This idea, originally due to Deutsch and Can- 
delas [ll|, was recently re-expressed by Graham, Jaffe 
and co-workers [2(| and by Barton [2l[. The implication 
is that in an experimental situation, a curved boundary 
might give rise to a large cutoff-dependent contribution 
to the physically measurable Casimir force [12] (if the 
boundary is not strictly rigid). This is consistent with 
Sen's observation that the coefficients of the cutoff- 
dependent terms contributing to the Casimir energy are 
geometrical objects such as the length of the circular 
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boundary. We also note that in the high-temperature 
limit the Casimir free energy is expressed in terms of a 
surface integral of a quadratic function of local curvature 
\'2-j which can be shown to have entirely geometrical na- 
ture. 

In order to better understand the geometrical meaning 
of the surface divergence we need to explicitly compute 
the Casimir energy due to a Dirichlet plane. This can be 
done by the method outlined in the previous Section. 

A. Surface energy of a plane in d dimensions 

Let us consider single Dirichlet plane perpendicular to 
the z-axis of the d-dimensional rectangular coordinate 
system. Since the space is uniform relative to transla- 
tions parallel to the boundary, both the fields u u (r) and 
f u (r±) (where is the position vector perpendicular 
to the z axis) can be expanded into a Fourier series; for 
example u u (r) = Y^q M uq( z ) ex P *Q r ±- The boundary- 
value problem ((4]) for the Fourier coefficients u u q(z) be- 
comes 

d 2 oj 2 

("To ~ q2 -)^q( z ) = °> u ^q(0) = /uq (12) 

dz z & 

where we assumed that the plane is located at z = 0. 
The particular solution (i.e. vanishing for / = 0) to (fT2"l) 
that decays as \z\ — >■ oo is 

^q(^) = /^qcxp (~vV + U) 2 / C 2 \z\} (13) 

Substituting the Fourier representations for u w (r) and 
fu(r±) along with the solution (fT3| in Eq.© we find 

S Z = W S ^ q 2 +U 2 /c 2 \f u q\ 2 (14) 
ui.q 

where A is the macroscopic (d — l)-dimensional area 
of the boundary. The action (fT4| has the required 
form Eq.©, with the geometrical coefficient Aq(w/c) = 
\j{2A\Jq 2 +oj 2 /c 2 ) that becomes small for large q. Ac- 
cording to Eq. (|13p the disturbance introduced by the 
fluctuating boundary condition is localized at the bound- 
ary, to within a length that is proportional to A itself. 
This is the source of the divergent coefficient of the sur- 
face tension. Substituting this A in Eq. (fTTj) and integrat- 
ing over ui we find 

4^ Jo " 2 2^ 2 

The sum in (|15j) can be recognized as the zero-point en- 
ergy of a harmonic field confined to a d — 1-dimensional 
space. By eliminating a range of the field modes, the 
Dirichlet boundary suppresses vacuum fluctuations in 
a region near it. This lowers the energy relative to 
the freely fluctuating vacuum, with the result that the 



Casimir energy due to the Dirichlet boundary (fl~5|) is neg- 
ative. The magnitude of the effect is a fraction of the 
d — 1-dimensional zero-point energy because the surface 
energy is dominated by the field modes highly localized 
at the boundary (without the cutoff this surface energy 
would be ultraviolet divergent). 

For macroscopic area A of the boundary the sum in 
(|15p can be transformed into an integral with the result 

e= jcAK^ r' q ^ dq (16) 

4 Jo 

where in taking the macroscopic limit we used the 
rule Y,q -> "4J d d - 1 q/(2Tr) d - 1 . The parameter K d is 
the surface area of a <i-dimensional unit sphere, Sd = 
27r d / 2 /r(d/2), divided by (2n) d : 

2-K d / 2 
Kd = {2n) d T{d/2) 

We note that for d > 1 the coefficient of surface tension is 
negative and given by 70 = £/ A ~ — fic(wo/c) d which is 
in accordance with the dimensional argument given ear- 
lier. For a smooth curved Dirichlet boundary, Eq. (|16[) 
represents the leading term of a geometric expansion to 
be discussed in more detail below. For d = 1 the sur- 
face energy (|16[) vanishes. This circumstance does not 
contradict our argument that introduction of the Dirich- 
let surface lowers the vacuum energy. It just means that 
in one dimension the leading finite-size correction to the 
bulk vacuum energy is negative, cutoff-independent, and 
of the order he divided by the system size. This correc- 
tion term is the difference between the sum (|15|) and the 
integral (1161) ; an explicit demonstration of this fact will 
be given below. 

B. Weyl expansion and the formally divergent part 
of the Casimir effect 

Although the surface energy is cutoff-dependent and 
determined by the transmission properties of the bound- 
ary, there is an intrinsic feature built into (1161) that can 
be revealed by rewriting the expression for the surface 
energy in a form that makes connection to its origin as a 
zero-point effect: 

/•oo Yicq 

£ = J —9area(q)dq (18) 

where 

9area(q) = -^AK d ^q d - 2 (19) 

can be called the areal density of states; it is a geomet- 
rical (and therefore generally relevant) quantity. This is 
an instance of a common feature, apparently first recog- 
nized by Dowker and Kennedy [lfj and by Deutsch and 
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Candelas and recently emphasized by Schaden [T3| 
that all cutoff-dependent contributions into the Casimir 
effect have a geometrical nature interpretable in terms of 
the DOS. Specifically, in the case of the scalar Casimir ef- 
fect considered in this paper, they all can be understood 
as having their origin in the asymptotic limit of the den- 
sity of eigenvalues of the Laplacian. The relationship to 
this classic problem of mathematical physics posed by 
Lorentz, tackled by Weyl and other researchers and 
made popular by Kac 2J| can be summarized as follows: 
The zero-point energy of the field confined to a region 
is the sum of zero-point energies of the field oscillators 



l hcq v 



Hcq 



G(q)dq 



where —q 2 are eigenvalues of the Laplacian: 
(A + q 2 )w = 



(20) 



(21) 



The discrete eigenvalue spectrum is determined by im- 
posing the Dirichlet boundary condition w(r) = at the 
boundary of the region. The function 



(22) 



is the exact DOS. It can be represented as a histogram 
with the distance between neighboring peaks of the order 
of the inverse of the system size. Quite generally [l2| the 
exact density of eigenvalues can be presented as a sum 
of the smooth (semi-classical) part g(q), commonly re- 
ferred to as the Weyl DOS, and an oscillatory remainder 
G(q)—g(q). The Weyl DOS g(q) can be understood as the 
result of averaging of the exact DOS (|2"2"|) over scales in 
q-space that exceed the distance between the neighboring 
peaks of G(q) [121] . Therefore g(q) contains information 
about the large-t? behavior of the exact DOS, and the 
cutoff-dependent part of the Casimir energy can be ex- 
plained as having its origin in the Weyl DOS. Likewise, 
the remainder G(q) — g(q) is responsible for the intrinsic 
part of the zero-point energy. It is then expected that 
in the macroscopic limit the zero-point energy (1201) nat- 
urally splits into cutoff-dependent and intrinsic pieces: 



£ = 



^9(q)dq 



g(q)]dq (23) 



In view of its ultraviolet (local) origin, the first cutoff- 
dependent term in ([2"B")) is additive, i.e. it can be pre- 
sented as a sum of volume, surface, line etc. boundary 
integrals. This is however not the case for the second 
(cutoff-independent) term where we have taken the limit 
of infinite cutoff frequency; requiring that the integral 
converges is a part of the definition of the Weyl DOS 

The smooth part of the exact DOS, g(q), can be repre- 
sented as a large-g expansion and each term of this Weyl 
expansion, similar to Eq. (|19p . can be interpreted geomet- 
rically. Indeed, for a d-dimensional volume V enclosed by 



a (d — l)-dimensional Dirichlet boundary of area A the 
Weyl expansion starts out as 

g(q) = VK d q d - x - l -AK d ^q d - 2 + ... (24) 

where the first term is the well-known consequence of 
the transformation rule ~ * ^ / d d q/(2ir) d . For an 
infinite vacuum only the first term is present; when this 
is the reference state the Casimir energy can be calcu- 
lated from Eq. (|23l) by removing the leading term from 
Eq. (f24|) . The second term, proportional to the area A 
of the boundary, only includes the effects of modes in- 
side the region and thus is half as large as (flU]) where the 
boundary has two sides. For d = 2, 3 the terms displayed 
in (jM]) are Weyl's original results [H, HJ]. For general d 
an equivalent expression for the number of states J g(q)dq 
has been given by Brownell [2~5j ] . 

Even though the Weyl DOS is a purely geometrical 
concept having little to do with physics, its relationship 
to the Casimir problem explains the sign of the surface 
term in Eq. (j24|) . 

A series of examples illustrating the general principles 
just outlined has been given by Schaden [TiJ. Below we 
consider two one-dimensional geometries that illustrate 
the expectation (|23|) . Additionally they contain the in- 
gredients needed in analysis of the scalar Casimir effect 
in two dimensions where we will find deviations from the 
rule 



C. One-dimensional interval 

Consider a harmonic field confined to a one- 
dimensional Dirichlet interval of length s. The eigenvalue 
spectrum is q n = irn/s, n = 1,2, ... and the exact DOS is 
given by 



G(q) = J2S(q 



-) 



(25) 



Since the peaks of this function are equidistant, the Weyl 
DOS can be readily obtained by noticing that every in- 
terval of length Aq = ir/s (that does not include the 
origin q = 0) contains exactly one eigenvalue. This im- 
plies that G(q) averaged over the interval of length Aq 
gives the Weyl DOS, g(q) = s/tt. This agrees with the 
general result (fMf and confirms that the surface (edge) 
term is indeed absent. 

The zero-point energy is given by 



TThC ^ , 

2s ^ 

n—1 



(26) 



In the macroscopic limit ui s/c 3> 1 that we are inter- 
ested in, the sum in (|26p can be computed with desired 
accuracy with the help of the Euler-Maclaurin summa- 
tion formula 12611 



— r°° ' 1 

]T 'F(n) a / F(x)dx - -F 

n=l •'° 



(0) 



1 

12 



F'(0) (27) 
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Applying this formula to the sum in Eq. (|26|) and employ- 
ing q = irn/s we find the expression 



£ = 



IT he 

~2s 



xdx 




hcq sdq 
~2 V 



irhc 
24s 



(28) 



that is in agreement with the general expectation (|23|) : 
the cutoff-dependent piece of the effect proportional to 
the system size s indeed originates from the Weyl DOS 
given by the d = 1 limit of (f2"4")l . g(q) = s/tt. The attrac- 
tive intrinsic part of the zero-point energy, — irhc/ (24s), 
is well-known [27j ; its sign was anticipated at the end of 
Section IIIA. 

The Casimir energy is known to be sensitive to the 
topology of the problem but this sensitivity is limited 
to the non-additive intrinsic part of the phenomenon be- 
cause it has non-local origin. As a simple illustration 
of this principle, let us consider a harmonic field con- 
fined to a one-dimensional interval of length s with pe- 
riodic boundary conditions. The eigenvalue spectrum is 
q n = 2im/ ' s, n = 0, ±1, ±2, ... and the zero-point energy 
is given by 



2irhc s-^ , 



2nhc 



xdx 

12 



hcq sdq 



it He 
~6s~ 



(29) 



where we applied the Euler-Maclaurin formula (|2T|) and 
employed q = 2irn/s. We observe that the first cutoff- 
dependent terms of (|28|) and (|29| are identical: both 
cases are described by the same Weyl DOS, and the dou- 
ble degeneracy of the \n\ > eigenvalues in periodic ge- 
ometry is compensated, compared to the Dirichlet case, 
by twice the distance between nearest g's. Only the last 
cutoff-independent terms of (f2"5)) and (f2T)]) sense the dif- 
ference in topology of the problems. The factor of four 
difference between the intrinsic parts of (|28|) and (l29l) is 
well-known [23] ■ The negative sign of the last term of 
(f2"9")l can be understood by realizing that the modes elim- 
inated by the periodic geometry do not contribute into 
the zero-point energy. 



D. Relationship to the zeta-regularized result 

At this point it is instructive to see how approach based 
on the concept of the Weyl DOS is related to regulariza- 
tion approaches that do not use physical cutoffs. To be 
specific, let us again consider the one-dimensional Dirich- 
let interval of length s. Then the ultraviolet divergence 
present in the cutoff-free energy sum (|26[) can be inter- 
preted in terms of the Ricmann £- function [5[ 



00 

£-< 

n=l 



(30) 



which is convergent for a > 1 and can be analytically 
continued to all complex a 7^ 1. Then Eq. ([2^|) is the 
value at a = — 1 of the expression 



irhc 

~2<r 



CM 



(31) 



where the superscript R stands for "regularized." Using 
the result £(— 1) = —1/12, this gives 



MR) _ 7rfic ^_n _ _^ hc 
t ~ 2s U j ~ 24s 



(32) 



which is the intrinsic part of (j28[) . The cutoff-dependent 
term of Eq. (|2"8"]) that is linear in s seems to have been dis- 
carded. However, we observe that it is possible to retrieve 
the cutoff-dependent term from the zeta-regularized re- 
sult as follows: the expression (|3"T1) has a pole at a = 1, 
which is between the range a > 1 where the sum (|3 1[) 
is convergent and the value a = — 1 where we wish to 
evaluate (|3~Tj) . We believe that in general there is a 
cutoff-dependent term associated with every such "by- 
passed pole" ; when the zeta-regularized expression is to 
be evaluated at a = — 1 and has a pole at a — m with 
residue A, the corresponding cutoff-dependent term has 
form J q m dq with a coefficient that is proportional to A. 
Then the zeta-function regularization method would give 
an elegant route to the determination of all of the terms 
in the Casimir energy. In fact Elizalde demonstrated that 
progress in this direction can be made if the zeta-function 
technique is supplemented with the Hadamard regular- 
ization method 1281 . 



E. Weyl energy and geometry of the boundary 

A natural consequence of the separability of the Weyl 
and intrinsic pieces of the Casimir energy is the conclu- 
sion that the Weyl energy £ (ujq) given by the first term 
of (f23| can be expressed in the form of a purely geomet- 
rical expression. When the media on both sides of the 
boundary have the same speed of light (so that locally 
the boundary cannot distinguish the inside from the out- 
side), the Weyl energy is given by a surface integral of 
an "even" combination of curvature invariants that does 
not depend on the sense of the local normal (the contri- 
butions from the "odd" terms cancel). For example, for 
a smooth boundary in three dimensions £(a>o) has the 
form [Till 



£(wo) 



ds(~f + 7i Q (Ci - C 2 ) 2 + 71&C1C2) (33) 



where 70 is the coefficient of surface tension mentioned 
earlier while 7i ,b are the curvature stiffness constants 
and Ci,2 are the principal curvatures. Since the boundary 
is made of real material, the shape coefficients 7„ have 
to be interpreted as contributions to the elastic constants 
of the boundary viewed as a flexible membrane. Assum- 
ing locality, the energy expression (I33[) can be written 
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down phenomenologically without referring to the Weyl 
problem. The latter however puts (|3"3")l on a solid foot- 
ing. Like the coefficient of surface tension 70 ~ Itiuq/c 2 , 
the curvature stiffness constants 7i a ,6 are dominated by 
highly localized field components that are suppressed by 
the boundary. Thus they can only depend on the mi- 
croscopic energy Hloq and the length scale c/wo; they are 
properties of the material of which the boundary is made, 
and are the same for all objects of whatever shape made 
from this material. Dimensional analysis then implies 
that Ji a ,b ~ fuJo We now see that the Weyl energy 



(|33[) is an expansion in powers of the cutoff frequency 
cjo with the leading uiq term and sub-leading luq term. 
In the macroscopic limit that we consider, there is no 
need to take into account powers of curvature invariants 
higher than those included in Eg. (1531 because the cor- 
responding curvature constants vanish in the luq — > 00 
limit. We also note that the geometrical expression (|3"3")l 
is applicable to an arbitrary harmonic theory and arbi- 
trary boundary conditions; however, the precise values of 
the stiffness constants 7's do depend on the nature of the 
field and type of boundary conditions. 

For simple geometries the expression for the Weyl en- 
ergy (|33| allows us to draw conclusions regarding whether 
the Casimir self-stress experienced by the boundary is in- 
dependent of the cutoff. Indeed, the Casimir energy for 
a spherical shell of radius a is expected to have the form 

m 



£. 



sphere 



47T70O 



„ he 

-4tt7i 6 + # — 
a 



(34) 



where the first two terms follow from (|33[) while the func- 
tional form of the last term is dictated by dimensional 
analysis; the numerical prefactor # is determined by the 
nature of the field involved and by the boundary condi- 
tions the shell imposes on the field. For the case of the 
electromagnetic field the coefficient of surface tension is 
known to be zero, 70 = [29j]. This means that the pres- 
sure on the boundary is independent of the cutoff and is 
correctly determined by regularization methods. At the 
same time it is clear that such a behavior is an excep- 
tion rather the rule. Indeed for the scalar Casimir effect 
with Dirichlet or Neumann boundary conditions the co- 
efficient of surface tension is non-zero, and the Casimir 
pressure is dominated by the first term of (|3~4"|) . 

For another example let us consider the case of an in- 
finite cylindrical shell of radius a. In the limit that the 
cylinder length L goes to infinity, its Casimir energy per 
unit length along the cylindrical axis is given by 



£ 



cylinder 



27TO70 



27T7i a 



#5 



(35) 



As in the spherical case (|3~4"j) . the cutoff-dependent terms 
are determined by Eg. (f3"3"]) while the functional form of 
the last piece is fixed by dimensional analysis. We now 
see that even in the case of the electromagnetic field, 
70 = 0, the pressure on the shell surface is dominated by 
the cutoff-dependent 2nji a /a contribution. 



F. Why are even space dimensions special? 

The analysis conducted up to this point assumed sep- 
arability of the Weyl and intrinsic pieces of the Casimir 
effect (see Eq. ([23|) ). In order to see that this assumption 
breaks down in even space dimensions, let us consider 
a sphere of radius a in d spatial dimensions. Assuming 
separability, the Weyl energy of the spherical shell can 
be estimated to have the form 



M 



M 



£(u ) ~ 51 7l 



E„ „d-l-2n 



(36) 



n=0 



n=0 



where in the first representation a~ 2n stands for the non- 
vanishing even power curvature terms and a rf_1 repre- 
sents the surface area of the boundary; numerical factors 
have been dropped for clarity. The rigidity constants j n 
have dimensionality of energy /length d ~ 1 ~ 2n , and they 
are determined by the microscopic physics, i. e. by thug 
and c/luq, the microscopic energy and length scales, re- 
spectively. This implies 



M 



z — ^ r 



71=0 



UJ 0- s d-2n a d-l-2r, 
C 



(37) 



We note that the functional form of the expression for the 
Weyl energy £(u>o) is more restrictive than what would 
be implied by a naive dimensional analysis that involves 
the energy scale fkoo and the two length scales c/luq and 
a. 

The number of formally divergent terms M + 1 in the 
expression for the Weyl energy £{wq) is fixed by the con- 
dition d-2M > 0. For odd d there are (d + l)/2 terms 
and the least divergent of these is linear in ojq. For even 
d there are 5 + 1 terms and the least divergent term is 
independent of the cutoff ujq. This however contradicts 
the expectation that the Weyl energy only contains the 
cutoff-dependent parts of the Casimir effect. The phe- 
nomenological resolution of this difficulty lies in allowing 
for a logarithmic cutoff dependence that replaces ujq in 
the expansion. Indeed if the Weyl DOS expansion con- 
tains a term proportional to 1/q 2 , then the first term of 
(|23|) is logarithmically divergent at low energies. Physi- 
cally the only low-energy cutoff in the problem is set by 
the inverse system size I/a. If this has to be invoked, 
then locality built into the geometrical expression for the 
Weyl energy is violated. This implies that when d is even, 
the Casimir energy has the usual cutoff-dependent (geo- 
metrical) and intrinsic (implied by dimensional analysis) 
contributions, and in addition a contribution of the form 



He ^ uiqcl 



(38) 



with a universal amplitude and logarithmically weak de- 
pendence on the cutoff frequency ujq. In the zeta- function 
regularization approach, this logarithmic divergence ap- 
pears as a pole exactly at a — — 1 , where the regularized 
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expression is to be evaluated. The detailed analysis of 
the two-dimensional case conducted below confirms this 
expectation, and in addition reveals the presence of an 
extra contribution into the Casimir effect that is both 
geometrical and independent of the cutoff. 



the geometry of the boundary, of the form 



CASIMIR EFFECT IN TWO SPATIAL 
DIMENSIONS 



IV. 



The two leading terms displayed in Eq. (|24|) are not 
sufficient to completely understand the cutoff dependent 
part of the Casimir effect in more than one dimension, 
and a more accurate expression is needed. An efficient 
indirect method of generating higher order terms of the 
Weyl expansion is via the Dirichlet series [H, H3] 

/>oo 

K(t)=Y / ^M-Qlt)= eM-1 2 t)G(q)dq (39) 

whose t — )• behavior probes the large-g piece of G(q), 
i.e. the Weyl DOS g(q). The latter, as implied by (1591 . 
can be extracted from the inverse Laplace transform of 
K(t->0) Bill. 

Below we present a detailed analysis of the two- 
dimensional geometry and determine both the cutoff- 
dependent and intrinsic parts of the Casimir effect. As a 
by-product we extract the Weyl DOS g(q) to desired ac- 
curacy. Our main result is that in the macroscopic limit 
the Casimir energy due to a smooth Dirichlet boundary 
r is given by 



£ = 



+ 



l/S 

hay 
256tt 



hcq 
hcq 



Adq 
'~2tT 
dq 



128irq 2 



C 2 (s)ds 



C 2 (s)ds+U n 



(40) 



where S is a macroscopic length scale specific to the prob- 
lem in question, s is the length along the boundary mea- 
sured relative to an arbitrary reference point, C(s) is the 
curvature at location s, and 7 = 0.57722 is Euler's con- 
stant. The Weyl DOS associated with (gDJ) is 



t \ A 1 
9{q) = -^~ 128V 



C 2 {s)ds. 



(41) 



The first term is the d = 2 case of the general result (fT9|) ; 
it is proportional to the total length of the boundary 
A = J r ds. The second term proportional to integrated 
square of local curvature agrees with that inferred from 
the t —> expansion of Eq.([39| due to Stewartson and 
Waechter HO]. 

The third term of (|40|) is cutoff- independent, even 
though it is geometry-dependent. This is a special feature 
of two dimensions. The first three terms of (1401) accumu- 
late a part of the Casimir energy which is dependent on 



■'geom 



hcq 



Adq\ 



2tt 



ficln^ 

C 

256?r 



C 2 (s)ds. 



(42) 

The integral leading to the second term of (|4"2"j) is evalu- 
ated with logarithmic accuracy and the outcome is com- 
bined with the third term of (|4U| . The presence of the 
infrared cutoff (set by finite system size) in the marginal 
second term of (|4"0"]) is necessary; it is another devia- 
tion from the expectation (|23|) that is due to the two- 
dimensionality of the problem. We note that in the first 
term of (|4"0"|) the lower integration limit is strictly zero; 
the associated higher order finite-size effects are accumu- 
lated in the last term of (l40l) . 

A similar expression has been previously given by Sen 
Q . However the sign of the first leading term of (1421) was 
determined incorrectly; additionally Sen's expression is 
inadvertently missing a factor of n in the denominator of 
the second term. While the first and second terms of (|40|) 
are dependent on the cutoff, the second term has only a 
very weak dependence. 

The first two terms of (|4"01) have their origin in the ad- 
ditive Weyl DOS (l4"Tj) while the third cutoff-independent 
term is additive as well. This is an aspect in which the 
last term of (j40|) denoted U na is different: it represents a 
non-additive intrinsic part of the Casimir effect which is 
sensitive to the topology of the problem. Its functional 
form is dictated by dimensional analysis: it is of the or- 
der he divided by a macroscopic length scale whose exact 
form is determined by the details of the problem. If there 
is only a single length scale a, then U na ~ Kc/a; in sev- 
eral cases below we will be able to explain both the sign 
and magnitude of U na that are not determined by the 
dimensional argument. The third and the fourth term of 
(|40p have the same order of magnitude; in view of their 
different physical origin we keep these terms separate. 

For a circle the effects of curvature and circumference 
are not distinguishable. We need a geometry where the 
two are independent, so that we can generalize the result 
to the case of an arbitrary boundary. To this end, we 
treat the wedge-arc geometry closely related to the sector 
geometry first considered by Kac [24j | . 



Wedge-circular arc geometry in two dimensions; 
Dirichlet case 



Let us consider a harmonic field theory (fTJ) defined in- 
side an infinite wedge of opening angle /3 whose edges 
enforce the Dirichlet boundary conditions. We are inter- 
ested in the change in the zero-point energy due to the 
introduction of a circular Dirichlet arc of radius a and 
length s = j3a as shown in Fig. 1. 

The semi-circular ( = tt) version of this problem has 
been studied earlier [3l|. Similar to the circular case, 
a divergence was found which was not removed by the 
zeta-function regularization method. 
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p 



a 



Figure 1: Wedge of opening angle f3 with superimposed arc 
of radius a in two dimensions. 



The wedge-arc configuration with an arbitrary angle j3 
has been considered in Ref. (32j with the aim of revealing 
the regularities in the contributions to the heat kernel co- 
efficients due to the failure of the boundary to be smooth. 
The local characteristics of the vacuum have been dis- 
cussed by Sakharian and collaborators for a scalar field 
with Dirichlet boundary condition in general space-time 
dimension [33| and for the electromagnetic field for d = 3 

M 

In polar coordinates p and ip the boundary- value prob- 
lem @ for this geometry becomes 



- — (P— ) 
p dp dp 



i d 2 

p 2 dip 2 



-7T ]U U = 



u u (p = a,tp) = f^{ip),u u {p 1 ip = 0(/3)) = (43) 

Seeking the solution in the form (p, ip) — 
12^=1 Run{p) s'm(iTnip / f3) and expanding the boundary 
field into a Fourier series f^{ip) = Y^=\ funSm(imip/l3) 
gives the boundary-value problem for the radial Fourier 
coefficients Runip) 



( d 2 



1 d 



\ dp 2 p dp ^ c 2 

(44) 

The particular solution to (|4*4l vanishing at p = 0, oo is 



Runip) 
Rum(p) 



I ml/fj {\uj\p/c) 

Jun T / 1 l / \ ' P ^ a 

I-nn/f}(Ma/c) 

Sum. „ /i I / \ ' p > a 

Kirn/piMa/c) 



(45) 



where I v {x) and K v (x) are modified Bessel functions [35j j . 
We note that as n increases, the functions (j45|) become 
increasingly more localized at the arc. 

If we take the simultaneous (3 — > 0, a — > oo limit but 
keep the length s = j3a fixed, the geometry of Fig. 1 turns 
into that of an infinite strip of width s whose sides are 
connected by a straight Dirichlet bridge. Indeed if we 
substitute p — a + z in Eq. (|4"4"l) . take the (3 — > 0, a —> oo, 



s = /3a — const limit and employ the fact that the role of 
the one-dimensional wave vector is played by q = nn/s, 
Eq.fH]) will reduce to Eq.([T2]). The discussion of Sec- 
tion IIIA then implies that in the limit neglecting the 
curvature of the arc the Casimir energy is given by neg- 
ative of the half of the zero-point energy of the harmonic 
field confined to a one-dimensional Dirichlet interval (see 
Eq. lfTS')) '). The zero-point energy of the field confined to 
one-dimensional Dirichlet interval was computed in Sec- 
tion IIIC (see Eq. (f2"8")l ). With that in mind we infer that 
neglecting the curvature the Casimir energy is given by 



hcq 



sdq 
~2~7 



nhc 
48s 



(46) 



This illustrates the effect of the geometrical size in gen- 
eral result (|40[) : the first formally divergent term has its 
origin in the Weyl DOS given by the first term of (|41l) . 
which is proportional to the arc length. The second in- 
trinsic term of (|46p has the expected Hc/s form and it is 
not additive. 

In order to account for the effects of curvature we sub- 
stitute the Fourier representations for u UJ {p 1 ip) and fuif) 
along with the solution (l4"5j) in Eq.©. Performing the 
angular integration we find 



Se - 2T^ 



2/„ 



*/p(\u\a/c)K vn/ p(\uj\a/c 



(47) 



where we employed the Wronskian relationship 

K v {z)l' u {z) - I u {z)K v {z) = 1/z [35]. The action 

(l47l) has the expected form ([6]) and thus the Casimir 
energy is given by the rule (fTTj) 



£ = 



he x -> 

2?ra ^ 

n=l 



dx In ( 2x1™ (x)Ke» {x) 



(48) 



In the f3 ^ I limit this expression can be analyzed with 
the help of the uniform asymptotic v 3> 1 expansion of 
Debye tUll]: 



2i/(l 



y 2 ) l ' 2 I v {vy)K v {uy) = 1 
6 5 



lis 2 X 1 + ', 



{l + y 2 ) 2 ' (1 + y 2 ) 3 



(49) 



If only the lowest order term of the Debye expansion is 
kept, the calculation of the Casimir energy (14"51) nearly 
mirrors that of Section IIIA; the outcome, in fact, repro- 
duces Eq.([46|). 

Writing the Casmir energy as £ — £W -f £( 2 ) _|_ ..^ 
using the next order term of the (3 <C 1 expansion (j49|) 
and integrating we find 



£(2) 



hes 



oo 1 

E - 



2567ra 2 L — ' n 

n—l 



(50) 



We observe that this term appears to be the length of 
the arc s multiplied by the square of the curvature 1/a. 
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Comparing it with the definition of the zcta-function (|30j) 
we see that without the physical cutoff the series in (|5TJ)) 
corresponds to the pole of the zeta-function. 

In the macroscopic limit that interests us, the sum in 
(f50|) can be transformed into an integral with the help of 
the definition of Euler's constant 7 [351 ] : 



00 

E - 

71=1 



dx 



This brings Eq. ([50)) into the form 



g{ 2) = f°° heq 

Jtt/s 2 



sdq 
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(51) 



(52) 



where we employed q = im/s. Combining Eas. (|46[) and 
(|52[) we arrive at the expression for the Casimir energy 
of the Dirichlet wedge of arc length s and small opening 
angle (3 



chlct 



7T j S 
HCJS 
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hcq 



sdq 
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1287TOV 
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(53) 



The first three terms are the geometrical part, and agrees 
with the curvature expansion (|40j) which is valid for an 
arbitrary curve. Indeed, since Eq. (|53[) is valid in the 
/3 = s/a <C 1 limit, we can take the arc length to be 
infinitesimally small, s — > ds. Then the first two terms 
of (|53[) can be viewed as having their origin in the differ- 
ential DOS 



, , , ds C 2 (s)ds 

d9iq) = ~^~ I2W 



(54) 



where C = 1/a is the curvature of the arc. Linearity in ds 
implies additivity. Then integrating over the boundary 
we arrive at Eq. ([4"Tj) . Similarly, the third term of (|53|) 
generalizes to the third term of (j40|) . 

The accuracy of the Debye expansion (14T)1) , implies that 
Eg. (1531 approximately captures the whole ^ (3 ^ 2tt 
range. 



B. Wedge-circular arc geometry in two dimensions: 
periodic boundary conditions 

The effect of topology can be investigated by compar- 
ing the geometry just studied with its periodic counter- 
part. Consider the wedge-circular arc geometry as in Fig. 
1 but assume that the wedge is bounded periodically, 
rather than by a Dirichlet edge, so that the condition 
u(p, if) = u{p, if + (3) applies. The opening angle will 
again be assumed to be arbitrary. For j3 — 2n the geom- 
etry is that of a plane with a Dirichlet circle; for general 



f3 we can envision the domain of the field as being the 
surface of a cone, again with a Dirichlet circle at distance 
a from the apex of the cone [32j . 

The Casimir energy of this configuration can be in- 
ferred from the solutions of the Dirichlet case. Indeed, 
the boundary-value problem we need to solve is posed 
by Eq.(|4"3l with the boundary field satisfying the con- 
dition of periodicity, fuj(f) = fu{(f + /?)■ This implies 
U u {p,(fi) = J2 < ^=- 00 R ^n(p)exp(2Trinif/l3). The radial 
function i? UIl (/?) satisfies the same Eq. (j4"4"|) with im/f3 
being replaced by 2im/ f3. The Casimir energy can then 
be inferred from Eq. (l4"51) as 

he 



2na 



dx \h(2xIo(x)Kq(x)) 



+ ^J2'J dx\n(2xI lfL (x)K 2f >(x)) (55) 

where the first term is due to the fluctuating angle- 
independent Fourier component of the boundary field / w 
while the sum is a contribution from the angle-dependent 
components. We note that for /3 — 2n Eq. ([55|) reduces to 
its circular counterpart [37[; this reference, however does 
not discuss the issue of the physical cutoff. 

The integral entering the first term of (I55|) was numeri- 
cally evaluated by Sen Q and more accurately by Milton 
and Ng [38[ with the result 

hc f°° , , n T , s„ , s 0.08808ftc , . 
£0 = 7^/ dx\n2xh(x )K (x = 56 

This is the intrinsic part of the Casimir effect unique to 
the periodic geometry; it naturally has no dependence on 
the opening angle of the wedge. The he/ a form of ([561 is 
dictated by dimensional analysis while the negative sign 
can be understood as due to elimination of the angle- 
independent modes by the boundary. 

The second term of Eq. (1551) can be analyzed analo- 
gously to that of Eq. (|48|) . Similar to the Dirichlet case 
studied in Section IVA, the results of the lowest order De- 
bye expansion can be inferred from analysis of periodic 
one-dimensional interval conducted in Section IIIC. 

The result of the calculation employing next order term 
of the Debye expansion is nearly identical to (|52"]) . The 
only difference is that the lower integration limit in the 
first term is replaced with 2ir/s; this is a consequence 
of the relation q = 2-Knjs and a sign of sensitivity to 
topology. Finally the counterpart of Eq. (|5"3"|) has the form 



periodic 



hcq 

~Y 

hcq 
~2 



sdq 
'~2n 



27r/s 

hcjs 

256-rra 2 



sdq 



1287raV, 
nhc 0.0880/ic 
12s 2na 



(57) 



This expression is again in agreement with general re- 
sults, Eqs. (|40p - (|4"2"|) . and the differences between the 
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Dirichlet (l53l) and periodic (|57|) Casimir energies are due 
to different topoiogy. The first source of the difference 
lies in different non-additive pieces of the effect U na (see 
Eq. (|4"0")) ). The second source, unique to two-dimensional 
geometry, is due to the second sub-leading term of the 
Weyl DOS (|4T|). 



C. Application: which geometry has lower energy - 
Dirichlet or periodic ? 

Here we apply our general results to the situation 
where the physics is dominated by the cutoff-independent 
part of the Casimir effect. Earlier we demonstrated that 
the Casimir energy due to a smooth Dirichlet curve is 
generally cutoff-dependent and dominated by the geo- 
metrical contribution (|42j) generated by the Weyl DOS 
(|41jl . However in the difference in energy between geome- 
tries that are only topologically distinct, only the cutoff- 
indcpcndcnt parts contribute. Indeed in the macroscopic 
limit uj s/c S> 1 the difference between (|57|) and (f53|) is 
given by 



AS = S 



periodic Srjirichlet 



Trhe 
16s 



0.08808ftc hcs In 2 



2na 



2567ra 2 
(58) 

where the last term originates from the second term of the 
Weyl DOS (l4"Tj) . It is straightforward to realize that for 
any relationship between a and s the periodic geometry 
has the larger energy, AS > 0. In order to see this we 
minimize (|58p with respect to a, substitute the outcome 
back in (f58| and observe that for any s we have AS > 0. 



cutoff-dependent contributions to the Casimer energy are 
determined by geometry. 

This fails in even dimensions, because the last relevant 
term of the Weyl DOS expansion is proportional to 1/q 2 . 
The corresponding contribution to the Casimir energy is 
only logarithmically dependent on the cutoff frequency. 
Then the separation of the cutoff-dependent and intrinsic 
parts of the Casimir energy is no longer clear cut. Here 
the only guaranteed way to arrive at physically measur- 
able answer is to use a cutoff method that emulates the 
transmission properties of the boundaries. Following this 
path we demonstrated that in the macroscopic limit the 
Casimir energy due to an arbitrary Dirichlet boundary 
in two-dimensions is given by Eq. (|40|) . Here the concept 
of the Weyl DOS again plays a prominent role but sep- 
aration of the intrinsic part of the effect requires careful 
analysis. This explains why the zeta-function regulariza- 
tion method does not work for circular [§3, HI] and semi- 
circular [3l| boundaries in two dimensions. Although we 
only analyzed the two-dimensional case, we expect that 
the mystery of divergent Casimir self-stress in general 
even space dimension Q is solved similarly. 

The approach to the Casimir effect employed in this 
work focused on the energy content of field modes elim- 
inated by the Dirichlet boundary. The applicability of 
this method is not limited to the scalar field theory ((T|), 
Dirichlet boundaries, or zero-temperature limit. One of 
the additional benefits of our method is the ability to 
provide the physical insight necessary to predict the sign 
of the Casimir self-stress. 



V. SUMMARY 
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